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ABSTRACT 

This  report  describes  the  multiple  precision  version  of  the  interval 
arithmetic  package  documented  in  MRC  Technical  Summary  Report  #1755.  The  mul- 
tiple precision  version,  based  on  the  FORTRAN  multiple  precision  arithmetic 
package  of  Brent,  is  extremely  portable.  It  is  assumed  that  the  reader  has 
access  to  TSR  #1755,  which  provides  the  basic  documentation  for  the  interval 
arithmetic  package;  the  current  report  addresses  only  those  aspects  of  the 
multiple  precision  version  which  are  not  covered  in  the  basic  documentation. 
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SIGNIFICANCE  AND  EXPLANATION 


Computer  calculations  are  inexact  due  to  rounding  error.  Interval 
arithmetic  is  a method  for  controlling  rounding  errors  in  such  a way  that 
an  answer,  say  x , is  expressed  as  a pair  of  numbers  (a,b)  between 
which  X must  lie  , a £ x _<  b . 

In  machine  interval  arithmetic  it  often  happens  that  interval  width 
b - a may  grow  rapidly  due  to  round-off  error.  Therefore,  in  some  computa- 
tions, single  precision  interval  arithmetic  may  not  be  adequate  to  produce 
satisfactory  results. 

This  paper  describes  a multiple  precision  version  of  the  interval  arith- 
metic package  documented  in  MRC  Technical  Summary  Report  #1755.  The  multiple 
precision  version  is  based  on  the  FORTRAN  multiple  precision  arithmetic  pack- 
age of  Brent,  which  is  described  in  ACM  Transactions  on  Mathematical  Software, 
volume  4 (1978),  pp.  57-70.  Because  of  the  portability  of  Brent's  package, 
the  multiple  precision  version  of  the  interval  arithmetic  package  is  also 
extremely  portable. 

The  "standard"  precision  of  the  multiple  precision  version  of  the  inter- 
val arithmetic  package  is  about  28  decimal  digits.  However,  the  user  may 
elect  to  use  a lower  precision  by  changing  certain  parameters  and  constants, 
or  may  alter  the  package  to  obtain  increased  precision.  This  report  explains 
the  basic  structure  of  the  multiple  precision  interval  package,  and  provides 
instructions  for  its  modification. 


responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
r.nnary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


1.  Introduction: 

I 

Because  interval  arithmetic  provides  rigorous  bounds  on  the  results  of  | 

computations,  it  is  necessary  when  implementing  the  interval  algorithms  on  a | 

computer  to  round  the  left  endpoint  to  the  nearest  machine-representable  num-  j 

ber  less  than  or  equal  to  the  true  result  of  the  calculation,  and  the  right 
endpoint  to  the  nearest  machine  number  greater  than  or  equal  to  the  true  re- 
sult. While  this  does  not  seem,  on  the  surface,  to  be  particularly  signifi- 
cant, in  long  computations  the  widths  of  intervals  may  grow  auite  large  due 
solely  to  round-off  error,  even  in  cases  where  such  inherent  problems  as 
dependency  are  not  present.  (The  classic  dependency  problem  arises  from  the 
inherent  inability  of  interval  mathematics  to  recognize  cases  where  the  same 
interval  appears  in  more  than  one  place  in  an  expression;  this  is  discussed 
elsewhere  (see,  e.g.,  [6]),  and  is  not  germane  to  the  present  discussion.) 

Thus,  in  some  applications,  single  precision  interval  arithmetic  may  not  be 
adequate. 

The  interyal  arithmetic  package  described  in  [7],  although  initially 
written  as  a single  precision  package  for  the  UNIVAC  1100  series  computers, 
was  designed  for  transportability  and  ease  of  modification.  When  the  multiple 
precision  arithmetic  package  of  Brent  ([1])  appeared,  it  seemed  logical  to 
proyide  a yersion  of  the  interyal  arithmetic  package  based  on  that  multiple 
precision  arithmetic  package.  Accordingly,  the  multiple  precision  arithmetic 
package  was  interfaced  with  the  AUGMENT  precompiler  ([3],  [4],  [5]),  upon 

which  the  original  yersion  of  the  interyal  package  depends;  this  interface  is 
described  in  [2].  The  interyal  arithmetic  package  was  then  adapted  to  use  the 
Brent  package  rather  than  machine  single  precision  format.  This  paper  de- 
scribes that  adaptation. 

In  Section  2,  we  describe  the  procedure  used  in  adapting  the  original 
package  to  produce  the  multiple  precision  yersion.  Section  3 describes  the 
procedure  for  bringing  the  multiple  precision  version  up  on  a host  computer 
other  than  the  UNIVAC  1100  series.  Section  4 contains  instructions  for  chang- 
ing the  precision  of  the  package,  both  statically  and  dynamically.  Section  s 
proyides  instructions  for  using  the  multiple  precision  yersion  of  the  interyal 
package.  The  Appendices  giye  the  AUGMENT  description  deck  for  the  multiple 
precision  yersion  of  the  interyal  package  (Appendix  1)  and  general  guidance 
for  setting  up  the  job  control  language  for  a production  run  of  a program 
which  uses  the  package  (Appendix  2). 

This  paper  is  intended  to  be  used  in  conjunction  with  the  basic  documen- 
tation for  the  interyal  arithmetic  package  [7].  Questions  concerning  the  op- 
erations and  functions  ayailable,  the  mathematical  basis  for  interyal  arithme- 
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tic,  the  archtecture  of  the  interval  package,  etc.  should  be  answered  in  that 
document.  Here,  we  address  only  those  aspects  which  are  unioue  to  the  multi- 
ple precision  version. 


2.  Adapting  the  interval  package: 

Five  separate  steps  were  required  to  adapt  the  original  interval  arithme- 
tic package  for  use  with  the  Brent  multiple  precision  arithmetic  oackage. 
These  were:  (a)  modifying  the  description  deck;  (b)  writing  the  BPA  primitive 
routines;  (c)  modifying  the  BPA  portion  of  the  package;  (d)  modifying  the  in- 
terval portion  of  the  package;  and  (e)  modifying  the  test  program.  These 
steps  are  discussed  in  greater  detail  below. 

Before  going  into  detail,  however,  a few  words  about  the  general  Philoso- 
phy of  the  adaptation  are  in  order.  Since  the  Brent  package  provides  not  only 
arithmetic,  but  all  appropriate  standard  mathematical  functions  as  well,  we 
wished  to  use  the  Brent  package  in  evaluating  the  standard  functions  for  in- 
terval arithmetic.  The  basic  interval  package  architecture,  however,  recuires 
that  routines  be  available  to  evaluate  the  standard  functions  in  higher  preci- 
sion than  that  of  the  interval  endpoints.  Moreover,  the  Brent  package  did  not 
provide  directed  roundings  in  the  basic  arithmetic  operations,  and  these  are 
required  by  the  interval  package.  Since  the  basic  arithmetic  operations  had 
to  be  rewritten  anyhow,  it  seemed  appropriate  to  use  the  Brent  package  in  the 
role  of  EXTENDED  arithmetic  (see  [?]),  making  the  interval  endpoint  arithmetic 
parallel  the  arithmetic  of  the  Brent  package,  but  with  fewer  significant  dig- 
its. Because  of  necessity  for  precise  input/output  conversions,  it  was  also 
regarded  as  desirable  to  restrict  the  exponent  range  of  the  interval  endpoints 
(BPA  numbers).  Accordingly,  we  take  the  BPA  numbers  to  have  three  fewer  sig- 
nificant digits  than  the  EXTENDED  numbers,  and  restrict  tne  exponent  range  to 
numbers  in  the  range  of  approximately  10*»-350  to  10»*350. 

These  remarks  having  been  made , we  are  now  in  a position  to  discuss  the 
details  of  the  modification. 

The  Description  Deck: 

The  modifications  to  the  description  deck  were  straightforward.  EXTENDED 
was  declared  as  type  MULTIPLE  rather  than  DOUBLE  PRECISION,  and,  since  this  is 
a nonstandard  data  type,  a copy  routine  had  to  be  named  in  the  declarations. 
BPA  was  then  declared  as  an  integer  array  of  length  9 (three  less  than  the 
declaration  for  MULTIPLE),  and  the  lines  which  named  unary  minus  and  copy  op- 
erators were  enabled  by  taking  them  out  of  the  comment  mode.  Field  functions 
allowing  access  to  the  sign,  exponent,  and  digits  of  BPA  numbers  were  added, 
and  conversion  functions  which  had  named  type  DOUBLE  PRECISION  were  changed  to 
name  type  MULTIPLE.  The  modifications  consisted  of  25  lines,  10  of  which  were 
editing  instructions. 


The  BPA  Primitives; 


These  consisted  of  the  four  arithmetic  operators  plus  conversion  from  EX- 
TENDED to  BPA  (with  directed  roundini?).  The  ADD,  SUBTRACT,  and  MULTIPLY  rou- 
tines were  patterned  after  Brent's  routines,  and  there  are  conseauently  sever- 
al supporting  routines  which  are  not  evident  (or  even  accessible)  to  the  user. 
The  rounding  routine  in  Brent's  package  was  altered  to  provide  directed 
roundings,  but  otherwise  is  quite  similar  to  Brent's.  The  DIVIDE  routine  had 
to  be  written  from  scratch,  since  Brent's  package  performed  division  by  takine 
the  reciprocal  of  the  divisor  and  multiplying;  this  is  faster,  but  does  not 
yield  the  precision  desired  for  BPA  arithmetic. 

The  only  other  primitive  required  for  this  implementation  was  a revised 
version  of  Brent's  MPINIT  — a routine  which  initializes  the  multiple  preci- 
sion package.  The  revised  version  also  initializes  the  BPA  package,  including 
setting  of  the  precision,  the  exponent  range,  and  calculation  of  the  constants 
needed  for  the  BPA  and  INTERVAL  portions  of  the  interval  package. 

Throughout  the  BPA  primitive  routines,  we  have  used  a COMMON  block  called 
3PAMCM,  which  corresponds  exactly  with  the  blank  COMMON  used  bv  Brent's  pack- 
age (except,  of  course,  that  BPAMCM  contains  the  corresponding  parameters  for 
the  BPA  numbers).  As  we  shall  see,  BPAMCM  must  be  altered  whenever  the  preci- 
sion of  the  package  is  to  be  increased. 

The  BPA  primitives  required  about  800  lines  of  FORTRAN  code;  however, 
these  primitives  should  not  need  to  be  modified  except  in  the  case  that  the 
precision  is  to  be  increased,  and  in  that  case  only  the  COMMON  declarations 
would  need  to  be  altered. 

Modifying  the  BPA  Portion  of  the  Package : 

The  modifications  to  the  BPA  portion  of  the  package  consisted  of  several 
parts;  The  conversion  arrays  had  to  be  expanded  to  accommodate  the  higher 
precison  and  larger  exponent  range;  the  constants  were  deleted  (since  they  are 
new  calculated  by  MPINIT) , and  various  routines  which  depend  on  EXTENDED  or 
BPA  formats  were  modified  — notably,  the  routines  which  pack  and  unpack  BPA 
numbers  and  convert  between  BPA  and  other  data  types.  In  all,  about  100 
lines,  approximately  one  quarter  of  which  were  editing  instructions,  were  re- 
quired uo  perform  the  modifications. 

Moc 1 fving  the  Interval  Portion  of  the  Package ; 

Aside  from  d'ileting  constant  definitions  (now  set  bv  MPINIT)  and  changing 
certain  routines  which  had  been  DOUBLE  PRECISION  FUNCTIONS  to  SUBROUTINES,  the 
majority  of  the  modifications  here  were  concentrated  in  the  INTRAP  routine, 
wi-cre  output  formats  and  lists  had  to  be  altered  to  conform  to  the  multiple 
P’^ecision  format. 

Two  subroutines  were  added  to  the  interval  portion  of  the  package.  The 
first  of  these  is  INTERS,  which  has  as  its  onlv  argument  an  integer  indicating 
thr  numoer  of  errors  that  are  allowed  before  execution  is  terminated  (unless 
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the  INTRAP  option  calls  for  unconditional  termiantion) . This  routine  is 
called  by  the  user’s  program,  if  it  is  used  at  all.  The  second  of  these  is 
INTWLK,  whose  task  it  is  to  print  a back-trace  of  subroutine  calls  in  the 
event  of  an  error.  This  routine  is  machine  specific,  but  must  have  two  argu- 
ments; the  first  being  an  identifier  for  the  calling  routine,  and  the  second 
an  integer  which  is  0 if  execution  is  not  to  be  terminated  and  nonzero  if  the 
error  counter  is  to  be  interrogated.  If  the  error  counter  is  interrogated,  it 
is  decremented  by  1 ; execution  is  terminated  when  the  error  counter  is  less 
than  or  equal  to  zero. 

In  all,  approximately  100  lines  were  required  for  the  modification;  about 
one  third  of  these  were  editing  instructions. 

Test  Program  Modification : 

About  140  lines,  approximately  50  of  which  were  editing  instructions, 
were  required  for  the  modification  of  the  test  program;  nearly  all  of  these 
were  concerned  with  I/O  formats.  Some  of  these  lines  will  need  to  be  revised 
if  the  static' precision  is  changed. 


3.  Bringing  the  Multiple  Precision  Interval  Package  up: 

The  first  step  in  implementing  the  multiple  precision  interval  package  on 
any  host  computer  is,  of  course,  to  obtain  a copy  of  Brent's  multiple  preci- 
sion arithmetic  package  with  the  AUGMENT  interface  (see  [2]  and  [1]).  Since 
this  package  was  not  developed  at  the  Mathematics  Research  Center,  it  is  not 
appropriate  for  MRC  to  distribute  it;  this  package  may  be  obtained  from 

Dr.  Richard  P.  Brent 
Computing  Research  Group 
The  Australian  National  University 
Canberra,  Australia 

Next,  one  must  have  a copy  of  the  AUGMENT  precompiler.  Unless  an  instal- 
lation already  has  AUGMENT,  the  request  for  the  multiple  precision  interval 
arithmetic  package  should  be  accompanied  by  a request  for  AUGMENT.  The  docu- 
mentation supplied  with  AUGMENT  will  include  the  necessary  instructions  for 
bringing  AUGMENT  up  on  the  new  host  system.  Both  AUGMENT  and  the  present 
package  may  be  obtained  from 

Mathematics  Research  Center 

The  University  of  Wisconsin  - Madison 

610  Walnut  St. 

Madison,  Wisconsin  53706 


With  AUGMENT  and  a copy  of  Brent's  multiple  precision  arithmetic  oaokaKe 
available,  the  next  step  is  to  bring  up  the  Brent  packaKe.  The  documentation 
for  that  package  includes  appropriate  guidance. 

The  multiple  precision  interval  arithmetic  package,  although  extremelv 
portable,  does  require  that  the  user  provide  three  primitive  routines:  the 
Hollerith  unpacking  routine  (a  similar  routine  must  be  provided  for  the  Brent 
package,  and  that  routine  may  be  used  for  the  interval  package  if  desired,  al- 
though it  must  be  modified  to  return  a buffer  length  of  133  if  the  unoacked 
string  would  exceed  132  characters);  the  error  handling  routine  INTERS  (for 
ERror  Set) , which  simply  records  the  number  of  nonfatal  errors  allowed  before 
execution  is  to  be  terminated;  and  the  routine  INTWLK  vrtiich,  in  the  UNI VAC 
1100  version,  prints  a trace  of  the  subroutine  calls  followed  by  a line  of  as- 
terisks. The  latter  routine  may  be  implemented  as  a routine  which  merely 
prints  a line  of  asterisks,  if  the  call  trace  is  not  easily  implemented  or  is 
not  wanted. 

Once  the  primitive  routines  have  been  provided,  it  is  necessary  only  to 
process  the  multiple  precision  interval  package  with  the  AUGMENT  precompiler 
and  then  compile  it  using  the  standard  FORTRAN  compiler,  placing  the  compiled 
modules  in  a library  file.  One  caveat:  we  found  it  necessary  to  increase  the 
working  storage  for  AUGMENT  in  order  to  provide  enough  work  space  to  handle 
the  descriptions  for  the  multiple  precision  package  and  the  multiple  precision 
interval  package;  5000  words  was  more  than  adecuate. 

The  test  program  should  also  be  AUGMENTed  and  compiled,  then  executed  to 
check  that  the  package  is  working  correctly. 

If  the  precision  of  Brent's  multiple  precision  arithmetic  oackage  is  to 
be  increased  over  that  defined  in  its  standard  AUGMENT  interface,  certain  mod- 
ifications will  also  be  necessary  in  the  multiple  precision  interval  package. 
See  Section  4 for  details. 


4.  Changing  the  precision  of  the  package: 
Static  precision; 


By  static  precision,  we  mean  the  maximum  allowable  precision  for  a multi- 
ple precision  interval  number.  This  precision  is  governed  by  the  number  of 
words  of  storage  allocated  to  each  BPA  variable;  this  is  9 in  the  standard  in- 
terface. Lower  precision  may  be  used  dynamically  (see  below),  although  if 
lower  precision  is  to  be  used  as  a standard  mode,  it  would  pay  to  revise  the 
package's  standard  precision.  Higher  precision  recuires  altering  the  dimen- 
sions of  certain  arrays.  Increasing  the  precison  may  be  necessary,  for  exam- 
ple, to  obtain  the  "standard  precision"  of  about  28  decimal  digits  on  short 
word-length  machines,  or  to  increase  precision  beyond  the  "standard"  oreci- 
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The  only  change  in  the  description  deck  is  in  Line  DSB00020,  where  the 
dimension  of  the  integer  array  for  a BPA  number  must  be  chaneed  from  9 to  n, 
where  n is  two  more  than  the  number  of  internal  base  digits  desired.  See  the 
documentation  for  Brent's  package  for  details.  The  only  thing  to  keep  in  mind 
here  is  that  n should  be  at  least  three  less  than  the  corresponding  value  for 
a multiple  precision  number. 

The  major  change  in  the  BPA  primitives  is  to  alter  the  length  of  the  R 
array  in  the  COMMON  block  BPAMCM.  An  array  of  length  4n  + 10  is  adecuate. 
This  array  is  referenced  in  the  following  lines  of  the  BPA  nrimitives: 

BPP00240 

BPP00810 

BPP01350 

BPP01980 

BPP02980 

BPP03990 

BPP04450 

BPP05220 

BPP06890 

If  an  exponent  range  larger  than  the  standard  one  is  desired,  line 
BPPO718O  must  be  changed  to  the  appropriate  exponent  (base  10).  In  this  case, 
do  not  neglect  to  make  the  appropriate  changes  in  the  conversion  arrays  in  the 
BPA  package  (see  below). 

The  COMMON  block  BPAMCM  also  appears  in  certain  rovitines  in  the  BPA  por- 
tion of  the  package.  Changes  must  be  made  in  the  following  lines: 

BPA01452 

BPA01724 

BPA05610 

BPA10580 

In  order  to  allow  working  space  for  the  input/output  conversion  of  an  ex- 
ponent greater  than  about  350  in  absolute  value,  the  lengths  of  the  conversion 
, arrays  ICX  and  ECX  must  be  increased,  and  the  position  of  the  radix  point  must 

be  altered.  See  [7]  for  details.  In  order  to  do  this,  chanees  must  be  made 
I in  the  following  lines: 

I BPA00430 

i BPA00480 

BPA00490 

I BPA01870 

I BPAO308O 

In  the  Interval  portion  of  the  package,  the  COMMON  block  BPAMCM  appears 
on  Line  INT1 1696.  Also,  in  Line  INTI 2920,  the  number  18  must  be  chanced  to 
2n,  and  in  Line  INTI 4290,  the  two  occurrences  of  916  must  be  chanced  to  nl6 
“ (if  n is  large  enough,  this  entire  format  statement  may  need  to  be  rewritten). 
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IT 


In  the  test  program,  BPAMCM  is  referenced  in  Line  TST1P400.  No  other 
changes  are  required  in  the  test  program,  although  many  of  the  formats  are  too 
short  even  to  display  the  full  significance  present  in  the  "standard"  oreci- 
sion . 


Dynamic  precision 


Dynamic,  or  run-time,  precision  may  be  altered  at  will,  within  the  bounds 
of  the  static  precision.  However,  there  are  several  things  to  consider: 


1)  The  precision  of  a BPA  number  must  be  no  greater  than  three  dieits 
less  than  the  precision  of  a multiple  precision  number. 


2)  If  precision  is  to  be  increased,  one  must  be  certain  that  the  digits 
beyond  those  already  used  in  each  BPA  number  (and,  by  extension,  each  interval 
endpoint)  are  zero. 


3)  If  precision  is  to  be  decreased,  special  care  must  be  taken  to  assure 
chat  the  new  intervals  contain  the  higher  precision  intervals.  This  is  true 
of  the  constants  in  the  package  as  well  as  those  intervals  that  may  have  been 
computed. 


In  this  event,  we  suggest  writing  a short  routine  which  will  round 
the  given  intervals  to  intervals  of  shorter  precision.  This  might  be  accom- 
plished as  follows: 


convert  the  left  endpoint  to  multiple  precision 


b)  decrease  the  precision  of  BPA  numbers  to  the  new  precision 


c)  set  OPTION  in  COMMON  block  BPACOM  to  the  value  of  RDL,  to  enable 
downward  directed  rounding. 


d)  use  BPACEB  to  round  the  left  endpoint 


set  unused  digits  of  the  left  endpoint  to  zero 


f)  repeat  for  the  right  endpoint,  using  the  value 
stop  corresponding  to  c). 


Once  these  precautions  have  been  taken,  change  the  value  of  T in  COMMON 
block  B ’AMCM  to  n-2,  i.e.,  the  number  of  digits  in  the  shortened  precision. 


5.  Using  the  multiple  precision  interval  package: 

The  only  unusual  feature  in  the  use  of  the  multiple  precision  interval 
package  is  that  the  user  must  place  the  statement 

INITIALIZE  MP 

before  any  statement  involving  multiple  precision  (EXTENDED),  BPA,  or  INTERVAL 
variables.  AUGMENT  will  translate  this  statement  to  a call  on  the  MPINIT  rou- 
tine, which  initializes  both  the  multiple  precision  package  and  the  interval 
package.  (Note  that  the  above  statement  will  cause  AUGMENT  to  declare  MP  as  a 
variable.)  Care  must  be  taken,  of  course,  to  insure  that  the  version  of  the 
MPINIT  routine  supplied  with  the  multiple  orecision  interval  packaee  is  the 
one  that  is  linked  with  the  program,  instead  of  the  version  suoolied  with  the 
multiple  precision  arithmetic  package  of  Brent. 

Otherwise,  the  use  of  the  multiple  precision  packase  is  exactly  the  same 
as  the  use  of.  the  single  precision  version;  see  [7]  for  details. 


6.  Conclusion: 

We  have  attempted  to  give  a concise  description  of  the  multinle  precision 
version  of  the  interval  arithmetic  package,  showing  how  it  was  derived  from 
the  single  precision  version,  how  it  is  implemented  on  another  host  system, 
and  how  it  may  be  used. 

This  document  is  intended  to  be  a supplement  to  [7],  and  that  report 
should  still  be  considered  the  primary  source  of  information  for  this  package. 
Naturally,  the  disclaimer  stated  in  that  document  applies  to  this  version  as 
well . 


Suggestions,  comments,  or  corrections  may  be  addressed  to  the  author. 
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APPENDIX  1 
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DESCRIPTION  DECK  FOR  MULTIPLE  PRECISION  INTERVAL  PACKAGE 


•DESCRIBE  EXTENDED 

DECLARE  MULTIPLE,  KIND  SAFE  SUBROUTINE,  PREFIX  EXT 
COMMENT  DUMMY  DESCRIPTION  SO  CONVERT  CAN  BE  USED 
SERVICE  COPY  (STR) 

•DESCRIBE  BPA 

DECURE  INTEGER(9),  KIND  SAFE  SUBROUTINE,  PREFIX  BPA 
OPERATOR  + (,  NULL  UNARY,  PRV,  $) 

OPERATOR  - (NEG,  UNARY,  PRV,  $) 

OPERATOR  + (ADD,  BINARY  1,  PRV,  $,  COMM),  • (MUL) , - ( SUB , , , , NONCOMM) , 
/ (DIV) 

OPERATOR  •»  (XBI,  BINARY  3,  PRV,  $,  INTEGER,  $) 

OPERATOR  .EQ.  (EQ,  BINARY  2,  PRV,  $,  LOGICAL),  .NE.  (NE),  .LT.  (LT), 
.LE.  (LE),  .GT.  (GT),  .GE.  (GE) 

FUNCTION  ABS  (ABS,  ($),  $) , SIN  (SIN),  COS  (COS),  TAN  (TAN),  ASIN  (ASN) 
ACOS  (ACS),  ATAN  (ATN),  L0G10  (LOG),  LOG  (LN),  LN  (), 

EXP  (EXP),  SINH  (SNH),  COSH  (CSH),  TANH  (TNH),  SORT  (SOT), 

CBRT  (CBT) 

FUNCTION  INT  (INT,  ($),  $) 

FUNCTION  MAX  (MAX,  ($,  $) , $),  MIM  (MIN) 

SERVICE  COPY(STR) 

FIELD  SGN  (rMPSIGA,  =MPSIGB,  ($),  INTEGER),  EXPON  (zMPEXPA,  rMPEXPB) , 
DIGIT  (=MPDGA,  =MPDGB,  ($,  INTEGER)) 

CONVERSION  CTB  (CEB,  MULTIPLE,  $,  DOWNWARD), 

CTB  (CRB,  REAL,  i,  UPWARD), 

CTB  (CIB,  INTEGER,  $,  UPWARD), 

CTE  (CBE,  $,  MULTIPLE,  UPWARD), 

CTR  (CBR,  $,  REAL,  DOWNWARD), 

CTI  (CBI,  $,  INTEGER,  DOWNWARD) 

•DESCRIBE  INTERVAL 

COMMENT  IN  THIS  DESCRIPTION  DECK,  'MULTIPLE*  IS  USED  AS  A 
SYNONYM  FOR  'EXTENDED'.  IF  EXTENDED  IS  A TYPE  OTHER  THAN  MULTIPLE, 
THE  APPROPRIATE  CHANGES  WILL  NEED  TO  BE  MADE  IN  THIS  DECK. 

DECLARE  BPA(2),  KIND  SAFE  SUBROUTINE,  PREFIX  INT 
OPERATOR  + (,NULL  UNARY,  PRV,  $) 

OPERATOR  - (NEG,  UNARY,  PRV,  $) 

OPERATOR  + (ADD,  BINARY1 , PRV,  $,  COMM),  • (MUL),  - (SUB, , , ,NONCOMM) , 

/ (DIV) 

OPERATOR  ••  (XXI,  BINARY  3,  PRV,  i,  INTEGER,  $),  ••  (XXB, , , ,BPA) , 

••  (XXE, , , ,DOUBLE  PRECISION),  ••  (XXX,,,, $) 

OPERATOR  .VEQ.  (VEQ,  BINARY  2,  .EQ.  , $,  LOGICAL,  COMM),  .VNE.  (VNE) , 
.SEQ.  (SEQ),  .SNE.  (SNE) 

OPERATOR  .VLT.  (VLT,  BINARY  2,  .LT.  , $,  LOGICAL),  .VLE.  (VLE), 

.VGT.  (VGT),  .VGE.  (VGE),  .SLT.  (SLT) , .SLE.  (SLE), 

.SGT.  (SGT),  .SGE.  (SGE) 

OPERATOR  .UNION.  (UNN,  BINARY  1,  .OR.,  $) , .INTSCT.  (SCT, , .AND.) 
OPERATOR  .SUBSET.  (SBS,  BINARY  2,  .EQ. , $,  LOGICAL),  .SPRSET.  (SPS) 
OPERATOR  .E.  (ELE,  BINARY  3,  .EQ. , BPA,  $,  LOGICAL) 


DSEOOniO 
DSE00020 
DSEn0030 
DSEoooio 
DSB00010 
DSBn0020 
DSB00030 
DSBOonyn 
DSBOOOSn 
DSBOOOqO 
DSB00100 
DSBOOI 10 
DSBOOlpn 
, DSBOOI 30 
DSBOOI UO 
DSBOOI so 
DSBOOI 60 
DSB00170 
DSBOOI 80 
DSB00230 
DSBOO234 
DSBOO236 
DSB00240 
DSB002S0 
DSB00260 
DSB00270 
DSP00280 
DSB00290 
DSX00010 
DSX00020 
DSXOOO3O 
DSX00040 
DSX00080 
DSXOOO9O 
DSX00100 
DSX001 10 
DSX00120 
DSXOOI3O 
DSX0O1UO 
DSX001SO 
DSX00160 
DSXOOI7O 
DSX00180 
DSXOOlpn 
DSX00200 
DSX00210 
DSX00220 
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FUNCTION  ABS  (ABS,  ($),  $)  , SIN(SIN),  COS  (COS),  TAN  (TAN),  ASIN  (ASN),  DSX00230 
ACOS  (ACS),  ATAN  (ATN),  L0G10  (LOG),  LN  (LN),  EXP  (EXP),  DSX0024n 

SINH  (SNH),  COSH  (CSH) , TANH  (TNH) , SORT  (SOT),  CBRT(CBT),  DSX00250 

LOG  (LN)  DSXn0260 

FUNCTION  ATAN2  (AT2,  ($,  $) , $)  DSX0n270 

FUNCTION  OK  (OK,  ($),  LOGICAL),  BAD  (BAD)  DSX00280 

FUNCTION  COMPOS  (CPS,  (BPA,  BPA) , $)  DSX00290 

FUNCTION  INT  (INT,  ($),  $)  DSX00300 

FUNCTION  DIST  (DST,  ($,  $),  BPA)  DSXOO3IO 

FUNCTION  MDPT  (MDB,  ($),  BPA),  HLGTH  (HLB),  LENGTH  (LGB) , LGTH  (LGB),  DSX00320 

MAG  (MAG),  MIG  (MIG),  PIVL  (PVL) , PIVU  (PVU) , SIZE  (SIZ)  DSX00330 

FUNCTION  SGN  (SGN,  ($),  INTEGER)  DSX00340 


FUNCTION  BOUND  (BND, 

(DOUBLE  PRECISION,  INTEGER),  i) 

DSXOO35O 

FIELD  SUP  (SUP,  SPL, 

($),  BPA),  INF  (INF,  INL) 

DSXOO36O 

SERVICE  COPY  (STR) 

DSXOO37O 

CONVERSION  CTX  (CIX, 

INTEGER,  $,  UPWARD), 

DSXOO38O 

CTX  (CBX, 

BPA,  $,  UPWARD), 

DSXOO39O 

CTX  (CRX, 

REAL,  $,  UPWARD), 

DSXOOUOO 

CTX  (CEX, 

MULTIPLE,  $,  UPWARD), 

DSX00410 

CTX  (ASG, 

HOLLERITH,  $,  UPWARD), 

DSX0042n 

CTI  (CXI, 

$,  INTEGER,  DOWNWARD), 

DSXOO43O 

CTE  (CXB, 

$,  BPA,  DOWNWARD), 

DSX00440 

CTR  (CXR, 

i,  REAL,  DOWNWARD), 

DSxno4Bn 

CTE  (CXE, 

$,  MULTIPLE,  DOWNWARD) 

DSX00460 
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i RUNSTREAM  FOR  MULTIPLE  PRECISION  INTERVAL  PACKAGE 

[INVOKE  AUGMENT] 

<description  deck  for  Brent's  packaRe>  j 

<description  deck  for  multiple  precision  interval  packaKe>  «.  I 

•BEGIN 

I •CONVERT  EXTENDED  - MULTIPLE 

<main  pro(5rain>  (does  not  apply  if  package  is  beine  comoiled) 

<nonexecutable  statements> 

INITIALIZE  MP 
<executable  statements> 

<subprogram(s)> 

•END  ; 

[COMPILE  PREPROCESSED  PROGRAM  MODULES]  ! 

<any  other  processing,  such  as  compilation  of  other  modules,  ! 

file  preparation,  etc.>  ■ 

[INVOKE  LINKAGE  EDITOR] 

<special  instructions  to  include  BLOCK  DATA  modules 
BPACOM  and  INTCCOM> 

I [ EXECUTE  PROGRAM] 

I <data>  (if  any) 

[END  OF  JOB] 

I 

t 

i 


♦ 


* <*• 


security  classification  of  this  page  (Wh««  0««  Bnfr^ 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


|2.  GOVT  accession  NO.I  3.  RECIPIENT'S  CATALOG  NUMBER 


REPORT  DOCUMENTATION  PAGE 


. REPORT  NUMBER 

1908  


4.  TITLE  faiHf  SuMltf*> 

THE  INTERVAL  ARITIMETIC  PACKKE  - MULT  1 1 
PRECISION  VERSION 


7.  AUTHORf*; 

J.  M.  Yohe 


S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Mathematics  Research  Center,  University  of  ^ 
610  Walnut  Street  Wisconsin 

Madison,  Wisconsin  53706 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U.  S.  Army  Research  Office 
P.O.  80x  12211 

Research  Triangle  Park,  North  Carolina  27709 


5.  TYPE  OF  REPORT  S PERIOD  COVERED 

Summary  Report  - no  specific 
reporting  period 


S.  PERFORMING  ORG.  REPORT  NUMBER 


a.  CONTRACT  OR  GRANT  NUMBER^*; 


DAAG29-7  5-C-0024  J 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  a WORK  UNIT  NUMBERS 


8 - Computer  Science 


12.  REPORT  DATE 

January  1979 


13.  NUMBER  OF  PAGES 
12 


4.  MONITORING  iGENCY  NAME  A AODRESSfff  from  Cmtrotlint  OHIem)  IS.  SECURITY  CLASS,  (ol  thit  rtport) 

UNCLASSIFIED 


tSa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (ot  thIt  Ktporl) 

Approved  for  public  release;  distribution  unlimited. 


t7.  DISTRIBUTION  STATEMENT  (oi  thm  mb^trset  in  Block  20,  It  dittcront  trom  Report) 


19.  key  words  (Continum  on  rororee  aide  II  noeoeeory  m%d  tdmntify  by  block  number) 

Interval  Analysis  Software  package 

Interval  Arithmetic  Precompiler  Interface 

' Multiple  Precision  Augment  Interface 

\ Extended  Precision 
' Portable  Software 


20.  ABSTI^ACT  (Continue  on  eereree  aide  If  neeeaeery  end  Identify  by  block  numbor^ 

This  report  describes  the  multiple  precision  version  of  the  interval  arith- 
metic package  documented  in  MRC  Technical  Summary  Report  #1755.  The  multiple 
precision  version,  based  on  the  FORTRAN  multiple  precision  arithmetic  package 
of  Brent,  is  extremely  portable.  It  is  assumed  that  the  reader  has  access  to 
TSR  #1755,  which  provides  the  basic  documentation  for  the  interval  arithmetic 
package;  the  current  report  addresses  only  those  aspects  of  the  multiple  preci- 
sion version  which  are  not  covered  in  the  basic  documentation. 


DO 

I JAN  73 


EDITION  OF  I NOV  61  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  RAGE  rWTiMi  Dmim  ErOrad) 


